Quantum entanglement in states generated by bilocal group algebras 



Alioscia Hamma, Radu Ionicioiu, and Paolo Zanardi 

Quantum Information Group, Institute for Scientific Interchange (ISI), Viale Settimio Severn 65, 1-10133 Torino, Italy 

Given a finite group G with a bilocal representation, we investigate the bipartite entanglement in the state 
constructed from the group algebra of G acting on a separable reference state. We find an upper bound for the 
von Neumann entropy for a bipartition (A, B) of a quantum system and conditions to saturate it. We show that 
these states can be interpreted as ground states of generic Hamiltonians or as the physical states in a quantum 
gauge theory and that under specific conditions their geometric entropy satisfies the entropic area law. If G is a 
group of spin flips acting on a set of qubits, these states are locally equivalent to 2-colorable (i.e., bipartite) graph 
states and they include GHZ, cluster states etc. Examples include an application to qudits and a calculation of 
the n-tangle for 2-colorable graph states. 

PACS numbers: 03.65.Ud, 03.67 .Mn, 05.50. +q 



I. INTRODUCTION 

Entanglement is certainly the most striking feature of quan- 
tum mechanics. Other than for its conceptual importance, en- 
tanglement has been in the last years one of the key concepts 
in quantum information, where it is the main resource re- 
quired in many protocols of quantum computation and quan- 
tum cryptography, e.g., quantum dense coding, Shor's algo- 
rithm and teleportation 1 1]. Entanglement is also an increas- 
ingly important concept in many topics of condensed mat- 
ter physics, like superconductivity and the fractional quantum 
Hall effect HI 01 Moreover, entanglement has recently 
been used, in the context of quantum phase transitions, as a 
novel tool to gain insights on the structure of the zero temper- 
ature phase diagram of interacting many -body systems |5|]. It 
also proves to play an important role in understanding certain 
aspects of quantum field theory [6] and spin systems ; 
these include 1 -dimensional lattice models for XY fl lldl and 
Heisenberg models • 

However, calculating and classifying entanglement for a 
general physical system is a daunting task. There is no 
known measure which completely characterizes the entangle- 
ment properties of an arbitrary system. If we restrict to bi- 
partite entanglement of a system in a pure state, this task be- 
comes easier. It has been proved that there is an essentially 
unique entanglement measure, namely the von Neumann en- 
tropy of the reduced density matrix of one of the two sub- 
systems, S = — Tr (pa log 2 pa) lil LLLIl - Although this looks 
conceptually simple, the calculation can be computationally 
intractable even for simple systems. 

On the other hand, in several quantum many-body sys- 
tems the states of physical interests (e.g., the ground state) 
are highly symmetric and these symmetries impose additional 
constraints which can simplify the calculation of entropy. This 
naturally leads to the topic of this article, namely exploring 
the entanglement entropy of a bipartite system using a group 
theoretical framework. 

Given a group G, possibly non-Abelian, with a bilocal ac- 
tion on a Hilbert space TL = Ha^'Hb, we introduce a class of 
states, the G-states, constructed from the group algebra of G 
acting on a (separable) reference state. G-states emerge natu- 
rally as the ground states of generic Hamiltonians. For these 



states we obtain an upper bound for the entanglement entropy, 
provided a separability conditions holds for the coefficients. 
A particular instance of this class, the G-homogeneous states, 
is constructed as an equal superposition of group elements 
acting on the same reference state |0). As expected, this ex- 
tra symmetry puts more constraints and is able to provide us 
an exact formula for the entropy. Remarkably, when A and 
B represent spatial complementary regions, e.g., in a lattice, 
we recover the entropic area law: the entanglement entropy 
for a bipartition (A, B) depends only on the degrees of free- 
dom localized on the boundary between A and B, and not on 
the bulk ones. This can be viewed as a another manifesta- 
tion of the Holographic Principle 1 12]. The entropic area law 
has been recovered in several physical systems H13L Il4l Hal . 
We examine several examples. If G is a group of spin flips 
acting on qubits, we show that the corresponding states are 
Calderbank-Shor-Steane (CSS) states flrjll and locally equiva- 
lent to a well-known class of stabilizer states, i.e., 2-colorable 
graph states. These states are important and appear in sev- 
eral physical contexts; examples include the ground state of 
the Kitaev model of topological quantum computation fnll . 
error correction codes (the CSS states) and the well known 
GHZ states. For arbitrary stabilizer states, entanglement has 
been studied also by Fattal et al. IU8I1 . This group theoretical 
framework proves to be fruitful also beyond the analysis of bi- 
partite entanglement for qubits. The first extension is to higher 
dimensional Hilbert spaces, i.e., qudits. The second one ap- 
plies to multipartite entanglement. Specifically, we show how 
to calculate the n-tangle for a G-homogeneous state, where G 
is a group of qubit spin-flips. 

The plan of the article is the following. The framework 
and the general bounds for the entropy in G-states and G- 
homogeneous states are given in Sectionln] These are a gen- 
eralization of our previous results IU91I20 1. Next we show that 
we can interpret the G-states as the physical states in a gauge 
theory (SectionllTH and we give a geometric interpretation of 
the entropic bound for G-homogeneous states and derive the 
area law. Section llVl is devoted to G-homogeneous states for 
a group of spin flips and we show that they are locally equiv- 
alent to 2-colorable graph states. We also discuss an example 
for qudits (Section flV C> and a calculation of multipartite en- 
tanglement (the n-tangle) for 2-colorable graph states (Section 
IvY We conclude in Section lvTI 
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II. ENTANGLEMENT IN G-STATES 

Consider (a unitary representation of) a group G acting on 
a Hilbert space TL. We assume that for a given repartition 
(A, B) of the Hilbert space, TL — TLa ® T~Lb, the action of 
G on TL is bi-local, i.e., any g G G is a bi-local operator 
.9 = 5a ® ,9b, where .9,4.5 are linear operators acting on the 
Hilbert space TLa.b, respectively. We also assume that it ex- 
ists a reference product state |0) = |fLi) <£> \0b) G W. We 
define the G-state as 

|* G > :=^a(. 9 ) 5 |0) (1) 

withE ff , fceG «( 5 )a(/i)(0|/i- 1 5 |0) = l. 

We point out that G-states appear naturally as the ground 
states of Hamiltonians belonging to the group algebra of G, 
H C C(G). Let D — {gi, ...,gk} be the set of generators of 
G with a local structure (i.e., each gi has a nontrivial action 
only on a set of neighbouring degrees of freedom). Then we 
can write any Hamiltonian in the form 

H = Y A (<?i)<?i (2) 

9i eD 

with X(g~ l ) = A(.g). Since the ground state of 
a generic Hamiltonian can be written as Vo) = 
lirn l g_, 00 Z~ 1 (/3)e~ l3H 0), then |^q) is a sum of all the ele- 
ments of the group with some coefficients A' (g) : Vo) — 
J2g£G nence it is a G-state. 

If all the coefficients are equal, we call the state a G- 
homogeneous state 

|G):=J\rVa£; s |0) (3) 

where the normalization factor TV := |G| ^ geG (0|.g|0) ^ Ois 
assumed to be nonzero and |G| is the order of G. It is obvious 
that the state |G) is stabilized by the group G, since h\G) = 
|G), V7i G G. Given a partition (A, B) of the Hilbert space 
TL we can write the state | G) as 

\G)=N- 1 I 2 Y j 9a\Qa)®9b\Vb) (4) 

where 3,4, b is the restriction of the operator g to the Hilbert 
space Ha,b- Consider now the two subgroups of G that act 
exclusively on the subsystems A and respectively B: 

G A ■= {geG \ g = g A ®^ B } (5) 
G B ■= {geG \ g = l A ®g B } (6) 

Let their order be d AlB '■= |Ga,b|- Since Ga and Gb are 
subgroups of the stabilizing group G it is true that g\G) = \G) 
also for every g G Ga.b- As these subgroups are normal, we 
can define the quotient group 

Gab := g^Tg-b (7) 



With this notation we have 

G = (J {{g A ®g B )h I g A ®t g G A , e G B ] (8) 

Proposition 1. Suppose we have a bipartition (A, B) of a 
Hilbert space, TL = Ha <8> Wg, and that the system is in a 
G-state |^g} = J2 g ec a (9)d\ty- We assume that for every 
g g G there is a representative /i, with [ft] g Gab, such that 
the coefficients satisfy the separability condition 

a(g) = a(g A ® 3bM = a A {gA)ctB{gB)(i{h) (9) 

where again su<g>l G G A , lgffs € Gb- Then the von Neu- 
mann entropy of the G-state corresponding to the bipartition 
(A, B) is bounded by 

S(|* G »<- Y \NAN B (3(h)\ 2 log 2 \N A N B (3(h)\ 2 

[h]<EG AB 

(10) 

where 

N 2 X := Y \ a *(9x)\ 2 , X = A,B (11) 

gx eGx 

Proof. From Eq.(|8j we have 

I*g) = Y a A (g A )a B {g B )(3(h){g A ®g B )h\0) 

g A ®g B £G A xG B 
[h]ec AB 

= (Qa®QbW g ) (12) 

where Q x ■= N x 1 Y /gxeGx a x {gx)gx, with X = A, B 
and 

|*' G ) - N A N B Y P( h ) h \°) 

[h]£G AB 

= N A N B Y PW\h A ) ^ \h B ) (13) 
[h]eG AB 

where \h x ) := h x \0x), X = A,B. Since Q A <8> Qb 
is a bi-local operator, the entanglement satisfies the bound 
SG^g)) < S(\t>' G )). The entropy of |*' G ) is maxi- 
mal when the set {l^x}}, X — A,B, is bi-orthogonal 
(Schmidt decomposition), in which case the entropy is 
-J2[h]eG AB \ N AN B l3(h)\ 2 log 2 \N A N B p(h)\\ hence this 
proves the bound. □ 

Observation. The separability condition seems rather 
strong. A simple example where this condition is satisfied 
is a(g) = x J (g)> i- e -, the coefficients correspond to ID- 
characters of the J-irrep of G. 

Corollary 1. If for any element g — gA ® 3b G G we have 

gx + lx => (0\gx\0) = (14) 
with X = A, B, then the entropy saturates its upper bound 

SQ*g)) = - Y WaNbPQi)? log 2 \N A N B p(h)\ 2 

[h]£G AB 

(15) 
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Proof. In order for the entropy 5(|^g)) to saturate the previ- 
ous bound we need to prove first, that the sets { | hx } , [h] G 
Gab} form a bi-orthogonal basis (hence S^l^g,)) saturates), 
and second, that \$>g) and I^g) have the same entropy. 

Consider now the scalar products {^[h^h'x |0) with 
h x ,h' x such that [h], [h'] G Gab and X = A,B. Then 
[h] ^ [ti] => h^h'x = (h^h^x # 1. (If = 1' then 

[fc- 1 /!'] = [1] = [h]- 1 ^'] [h] = [/i']). FromEq. d it 
follows that the sets {\h x ) \ [h] G Gab}, with X = A, B, are 
a bi-orthogonal basis. 

For the second part it is enough to show that the set 
{[hx) ■— Qx\hx), X = A, B} is bi-orthogonal. Consider 

(hx\h' x ) = (hx\Q x Qx\h'x) 

= N x 2 a x {gx)ax{g' x ) 

gx,g' x eG x 

x (Olh^g^g'xh'xlO) (16) 

By hypothesis i ll-It the last scalar product is zero unless 

h x 1 9~x9x h 'x = hence h x = (g^d'xWx and thus 
[hx] — [h'x] since they differ by an element of Gab- As we 
can choose the same representative for a given equivalence 
class, h' x = hx and hence g' x = gx- Taking into account 
the normalization (II 1> . we obtain (hx\h' x ) = 5j lx ^ ■ Thus 
\^g) and YS>' G ) have the same entropy and the thesis follows 
immediately. □ 
Corollary 2. For G-homogeneous states the entropy is 
bounded by 

S<log 2 |G AB | (17) 

and the bound is saturated if the condition (1 14-i holds. 

Corollary 3. If the group G is a direct product G = Ga x Gb, 
then S = 0. 

Proof. In this case we have Mg G G, g — gA <8 gs with 
gA <8 1b G G^ and 1U G Gb ■ Then we can write the 
G-state as 

|*g) = J! "a3a|0a)» X a s5s|0s) 

=: |XA>®|XB> (18) 

As this is a product state with respect to the (A, B) partition, 
it is obvious that its entanglement is zero. □ 

III. G-STATES, GAUGE THEORIES AND THE ENTROPIC 
AREA LAW 

What is a G-homogeneous state? And for which states 
the bound in Corollary 2 is saturated, i.e., S — log 2 |Gab|? 
The following construction can give us some insight. We 
assume that Ti has a given tensor product structure TL ~ 
C dl ® ... ® C dL . Let |0di) be a reference vector in C di and 
|0) the product state |0) := |0 dl ) ® ... ® \0 dL ). Now let us 
construct one possible set C := {Gi, Cl} of linear op- 
erators Gj G £(C dl ), Vt = In the following we 



take G to be a finite group of linear operators in TL gener- 
ated by a suitable choice of tensor products of elements in 
C : G — {Ai, ...,A n ) , where every 

A k = d (19) 

ie.Jk 

is a local operator on the local Hilbert space TLk = ®i e j k C di 
and Jk is a set of indices. We call the Ak star operators 
for reasons that will be clear in the following. A bipartition 
(A, B) is given by the choice of indices I a, Ib belonging to 
the subsystems A and B respectively, such that I a U Ib = 
{1, L} and TL — Ha ® H B , with 7i A = <S> ie i A cd '> and 
similarly for TLb- Notice that automatically all the operators 
in G have a bilocal action on every partition (A, B) of the sys- 
tem. If two sets of indices Jk , Jk' have one index in common, 
Jk n Jk' = {i}, we say that the corresponding star operators 
overlap on the local Hilbert space labeled by i. 

Let us now consider n sets of indices J\ , . . . , J„ defining n 
star operators such that: 

Jkc\Jj r , k = l,...,n (20) 

r^k 

VJ fe! J k ' : J k D J k ' ={i}or0 (21) 

We can see a graph (or a lattice) emerging from these n star 
operators. The vertices of the graph correspond to the star op- 
erators, so to n star operators we associate n vertices. The 
edges of the graph are the local Hilbert spaces C di on which 
the star operators overlap. This justifies their name: a star 
operator Ak is the product of the operators associated to the 
edges extruding from a vertex k. Star operators obviously 
have a local structure. 

With this construction, we obtain a geometric interpreta- 
tion of the entanglement in G-states. First of all, the en- 
tanglement with respect to a partition (A, B) is zero if we 

( x) 

can split the set {Ak} in two subsets {A^- } such that 

= ® ie j KC j x G, X = A,B. So the A k 's in one 
subset do not overlap with the ones in the other. The two sub- 
sets generate together the whole G and since each generates 
Gx , X — A,B, respectively, in this case G is a direct product 
G = Ga x Gb ; hence the entanglement is zero from Corollary 
3. Geometrically, this means that A and B are not connected 
in the graph. If the graph is connected then there is no parti- 
tion for which the entanglement is zero. On the other hand, if 
no A^s overlap, the graph is made of all isolated points and 
the G-state is completely disentangled. If we choose the Ak's 
such that we form a lattice, then there is no way to partition 
the system such that G = Ga x Gb and the entanglement is 
always different from zero under any possible partition. 

We can view the Ak's generating the group G as local gauge 
symmetries. The vertices of the graph represent points in 
space. Being G Abelian means that the symmetries act in- 
dependently at every point. In quantum mechanics we can 
construct a gauge theory by projecting a Hilbert space TL to 
a smaller Hilbert space TL p hys of the physical states. This is 
done by requiring that the physical states are annihilated by 
some operators; for example, the physical states \(f> p hy S ) in 
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quantum electromagnetism are the ones annihilated (at every 
space-time point) by d^A^phys) — 0. The Hilbert space H 
is the total Hilbert space. The physical states are the states 
annihilated by the operators Ak — 1 ll2lll : 

H ph ys = {0 E H\ {A k - l)il> = 0} (22) 

The G-homogeneous state is obviously a physical state. 

How can we characterize all the physical states? Let us 
find the algebra of the linear operators acting on the phys- 
ical Hilbert space C(H p h ys )- What does this algebra look 
like? Consider the subgroup W C C(H) of linear operators 
commuting element-wise with G and let T{W) be its asso- 
ciative algebra. Then consider the ideal I generated by the 
Ak — 1. The algebra acting on the physical Hilbert space is 
then C(TL p hys) = J-(W)/X. Then it is immediate to see that 

H ph ys = T{W)/1\G) = span{W|G}} (23) 

In other words, the orbit of W through the G-homogeneous 
state |G) is an orthonormal basis in Hphys- In general a phys- 
ical state is not an equal superposition of a group of linear 
operators, hence we cannot apply Proposition 1 to any quan- 
tum physical state of a gauge theory. Nevertheless, consider 
this particular situation: the system is partitioned such that ev- 
ery element of VV is gauge equivalent to some operator acting 
exclusively on A or B. Then, in the hypothesis of G and W 
belonging to the subspace B C C(H) of linear operators such 
that (0|7|0) = for all 7 G B, the following proposition holds 
for every physical state in Hphys'- 

Proposition 2. Let (A, B) be a partition of H such that W 
acts (modulo a gauge transformation g G G) exclusively on 
either A or B, say w ~ 1a <8> wb, Vu> G W. Then the 
entanglement entropy of every physical state \ip) G Hphys is 
equal to S = log 2 |Gab|- 

Proof. The physical states in the physical Hilbert space H p h ys 
can be written as 

Hphys) = C ( W ) W \ G ) (24) 
mew 

with the normalization condition 2~2 w ew \ c ( w )\ 2 = The 
density matrix of this pure state is 

pPhys = c(w')c(w)w' pw- 1 

= c(w')c(w~ i )w l pw (25) 

w,w' £W 

where p is the density matrix of the G-homogeneous state, 
p := |G) (G|. Tracing out the B degrees of freedom we obtain 
the reduced density matrix 

pPhys = cK)c(^ _1 )Tr B (w' 'pw) (26) 

By hypothesis for every element w G W there is a gauge 
transformation g G G such that w := gw = 1^ ® {gw)s- 
Denote g, g' such gauge transformations for the operators 



W, w' in Ea. J26l l. Then obviously wpw' = wgpg'w' = 
wpw', since [w, g] = and = p. The particu- 
lar form of w implies that Tr b (wpw 1 ) = Tr siw'wp) — 
unless w' = w^ 1 . This follows from Tr s{wp) — 
E S9 ' gg 9a\0a)(0a\9a9a 1 (0i3|ffsw|°B> and for w ^ lithe 
scalar product (0B|<7^t()|0B} is always zero. The reduced den- 
sity matrix becomes 

pj^ys _ c(w' ) c(w~ 1 )Tr b(w pw) 

w,w' £W 

= ^2 c(w')c(w' 1 )S li , {b -iTr B p 

= PA (27) 

Then the von Neumann entropy for every physical state is 
equal to the entropy of the G-homogeneous state p. □ 
Entropic area law. The entropy of the physical states in H p h ys 
satisfies the area law. Consider a partition (A, B) of H p h ys 
constructed by taking as subsystem A all the degrees of free- 
dom (i.e., the local Hilbert spaces C di corresponding to the 
edges of the graph) inside or intersected by a closed sur- 
face S. The group Gab will be generated by all the stars 
Ak based on sites outside the surface that puncture the sur- 
face E. Let tiab the number of such stars. Then the en- 
tropy is S = log 2 I Gab I = log 2 /(n^), where the function 
/(riAs) := I Gab I gives the order of the group as a func- 
tion of the number of its generators. As a measure of the sur- 
face E we can choose the number of punctures a. In general, 
o > n AB because a star can puncture the surface in more than 
one point. It immediately follows that 

S < log 2 f(a) (28) 

This bound can be saturated for some geometries, e.g., in a 
cubic lattice if we choose E to be convex 12011 . 



IV. G-HOMOGENEOUS STATES FOR SPIN FLIPS 
A. Qubits 

In our previous work II 9l 12011 we investigated examples of 
G-homogeneous states where G was a group of spin flips and 
hence the states were also stabilizer states; in this case G is 
Abelian and moreover Vg G G, g 2 = 1. We will clarify later 
the relationship between this class of G-homogeneous states 
and the well-known graph states. More exactly, we will prove 
that G-homogeneous states corresponding to a group of spin 
flips are locally equivalent to 2-colorable (i.e., bipartite) graph 
states. This representation of 2-colorable graph states as an 
equal superposition of elements of G acting on a reference 
state proves to be very effective in calculating entanglement 
and general correlation functions. 

In order to make the connection between G-homogeneous 
states and stabilizer states we start with some general con- 
siderations about the stabilizer formalism. Consider a sys- 
tem of n spins- 1/2 with Hilbert space H = Hf n (Hi = 
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span{|0), |1}} is the Hilbert space of a single spin). As be- 
fore, we choose the reference state |0) := |0)i ® • • • <8> |0}„, 
corresponding to the basis vector with all spins up. We de- 
note the Pauli matrices by Xi, Yi and Zi (the subscript rep- 
resents the qubit on which they act). There are two nec- 
essary and sufficient conditions for a group of Pauli opera- 
tors S = (si, . . . , s n ) to be a stabilizer group Jj], l26ll : (i) 
sf = 1, Si — J, Vi; and (ii) [si,Sj] = 0, Vi,j. From 
the general theory we know that any element of a stabilizer 
group can be written as ±X(a)Z(b), where a, b G ~%- 2 
are binary vectors and X(a) := Y[ k =iX k k ', analogously, 

Z(b) = nLi z k k l 22 l- We can define the nx 2n gener- 
ator matrix as follows 



ai 



b„, 



(29) 



where Si = ±X(a.i)Z(hi) and a^b^ G ~2- 2 are the binary 
vectors corresponding to the generators (the sign is omitted in 
the definition of the generator matrix). Hence the left (right) 
part of the generator matrix contains the X(Z)-type gener- 
ators, respectively. We also define the inner product of two 
binary vectors as 



(a, b) := aibi mod 2 

i=l 

The two conditions satisfied by the generators become: 

(i) (a,,b l ) = 

(m) (ai,bj) + (&j,bi) — mod 2, Vi, j 



(30) 



(31) 
(32) 

of 



An n-qubit stabilizer state is stabilized by a group 
Pauli operators having n generators. 

For example, the Bell states |$±) = 2~ 1 / 2 (|00) ± |11)) 
and |* ± ) = 2" 1 / 2 (|01) ± |10)) are stabilized by the groups 
(±XiX 2l ZiZ 2 ) and (±X 1 X 2l -Z\Z 2 ), respectively and the 



generator matrix is 



11 
00 



00 

11 



Let C\ — {H,P) be the 1 -qubit Clifford group, mapping 
products of Pauli matrices into products of Pauli matrices; 
H = 2- 1 / 2 (X + Z) is the Hadamard and P = diag(l, i). 

A natural question to ask is when two stabilizer states are 
locally equivalent and hence they have the same entropy S. A 
restricted criterion for local Clifford equivalence only is the 
following: 

Proposition 3. Let G and G' be two stabilizer groups of n 
qubits and let G be generated by G = (gi, <?„}. Then the 
states stabilized by G and G' are locally Clifford equivalent 
iff there is a local unitary U G Cf n such that the set hi = 
UgiU\ i = 1, n generates G 1 = (hi, h n ). 
Proof. Denote by pa and pa> the states stabilized by G and 
G', respectively. Since G is generated by (gi, g n ), we can 
write 



PG 



= \G)(G\=2- n Y[(l + gi ) 



(33) 



The two stabilizer states are locally Clifford equivalent iff 
there is a local unitary U € Cf 11 such that p G , = U p G W = 
2~"IlLi( 1 + Ug t W); define hi := U gi U\ i = l,...,n. 
Then hf — 11 and [hi, hj] = 0, Vi, j, since the generators gi 
satisfy the same relations. Then hi can be chosen as a set of 
generators for G' . □ 

Turning to the previous example, it is easy to see that 
|$~) = Zi|$ + ) and hence their generators are related by 
Zi(XiX 2 )Zi = —X\X 2 , Zi(ZiZ 2 )Zi = Z\Z 2 . 

Let G C {Xi, . . . ,X n ) be a group of spin flips: V<? £ 
G,g = ^(a), where a 6 W. 2 is a binary vector and X(a) = 
Yl"—i Xi*- Consider the generators of G = (gi, . . . , <j%) = 
(X(ai),...,X(a k )); \G\ = 2 k . Then we can write the 
group G in terms of an (additive) group of binary vectors 
A = (ai, . . . ,a k ) C Z£, and therefore G = X(A). 

For a group G of spin flips, the G-homogeneous state is: 



|G) = |G|- 1 / 2 ^X(a)|0) 



(34) 



aeA 



As mentioned before, G leaves invariant the state |G), since 
g\G) = |G),Vg G G. We now construct the stabilizer of 
\G),S\ G) = (gi,..,g k ,s k+ i,..,s n )\ obviously it should have 
n generators. The first k are the pure spin flips generating G, 
gi = X(ai), 1 < i < k; the other generators can be written as 
Sj = ±X(a j )Z(b j ) with k + 1 < j < n. Since Sj \G) = |G) 
and [sj,g] = 0, Vg, we have 



1 = <G| Sj |G)=±|G|- 1 (0\gX{a 3 )Z{h 3 )g'\0) 



g,g'eG 



±J2{0\hX(a j )\0)=± XG (X(a j )) 



(35) 



heG 



where the characteristic function is defined as xc{h) = 1 
if ft G G and otherwise. This implies that ±X(aj) G 
G, hence we can choose all the generators Sj as pure 
phase flips, Sj — Z(hj). Therefore the stabilizer 6>|g) = 
(.X-(ai), ..,X{a k ), Z(b k+1 ), .., Z(b n )) = X(A) ■ Z{B) is a 
direct product of pure spin flips and pure phase flips, with 
A = (ai, .., a k ) and B = (bfc+i, ■■, b„). Commutation of all 
generators implies (a^bj) = 0, Mi, j, hence B = A^. Thus 
the whole information about the state is contained only in the 
AT-type generators (or, equivalently, only in the Z-type ones): 



S lG) = X(A) ■ Z(A^) 



(36) 



and so the group G = X(A) describes fully the stabilizer state 
|G). A stabilizer state having only X- and Z-type generators 
is also known as a Calderbank-Shor-Steane (CSS) state. 

Another way of seeing this is the following. The stabilizer 
of the vacuum is <S|o) := (Z\ , . . . ,Z n ) and hence: 



2-" h = 2 ~" z ( b ) 



(37) 



hes 



As any group of pure phase flips Z(B) C 5|o), these sym- 
metries are already included in the way we constructed the 
G-homogeneous state P4| |. since the reference state is exactly 
|0>. 
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B. Relationship to graph states 

Graph states are a class of multiparticle entangled states 
which include cluster states, GHZ states etc and have been 
extensively studied recently B231I24H25II. Given a graph Q, the 
associated graph state is stabilized by the following n genera- 
tors: 



fj, 



=x t n 



i = l. 



(38) 



where the product is taken over all the nearest-neighbor ver- 
tices of vertex i. It is easy to see that the n x 2n matrix de- 
scribing the generators is the following: 



(1„|M) 



(39) 



hence the X-part of the generators is the unit matrix and the 
Z-part the adjacency matrix M of the graph. Let us analyze 
first the relationship between stabilizer states and graph states. 
Clearly, graph states are a particular class of stabilizer states. 
According to Gottesman ([26], Chapter 4) we can put the sta- 
bilizer matrix into a standard form by performing Gaussian 
elimination. For stabilizer states (set k = in Eq.(4.3) of 
1 26]) we have the standard form 



Jl, 



B r 



Cr.n — r 



(40) 



where r is the rank of the X part of the generators matrix (the 
subscripts denote the size of the matrices). If r = n, then this 
is equal to 



(41) 



with both n x n-matrices. Now, is B the adjacency matrix of 
a graph? We need to prove that Bu = and By = Bji. This 
follows immediately from the properties of stabilizer gener- 
ators. For any generators gij = X(a.ij)Z(hij), we have 
(aj,bj) = and (a^b,-) + (aj,bj) = mod 2, Vi,j',from 
which follows that B is indeed an adjacency matrix. There- 
fore, if the rank of the X part of the generator matrix is max- 
imal, r — n, the stabilizer state in question is a graph state. 
What if r < nl It has been shown that any stabilizer state 
is equivalent to a graph state under local Clifford operation 
llz^l . In conclusion, an arbitrary stabilizer state is either a 
graph state (if r — n), or is locally equivalent to a graph state 
(r < n). 

One problem with graph states is that there is no 1-to-l cor- 
respondence between graphs and locally inequivalent states. 
Thus, a GHZ state can be described either by a star graph or 
by a fully connected graph |28|. Moreover, some graph states 
have a more compact description as G-homogeneous states. 
For example, the stabilizer of the GHZ state is 



Sghz = {X n , Z\Z%, ... , Z±Z n ) 



(42) 



with the notation X® n := Yl7=i Xi- Described as a graph 
state, the GHZ state is (locally equivalent) to a star graph with 
a stabilizer group generated by 



S' star — (Xl J^J Zi, X2Z1, X n Zi) 



(43) 



It is immediate to see that the two set of generators are related 
by g' k = HgfeH, k = 1, ...,n, where the local operator is 
H := rifc=2 Hk (Hk is a Hadamard on the k qubit). Then the 
density matrix of the star graph state stabilized by S' s tar is 



p star = 2- n (l+X 1 ]JZ i )l[(l 

i=2 



■XjZx) (44) 



i=2 



The above form is more complicated than the description of 
\GHZ) in our formalism, since the group G of spin flips is 
just G = (X® n ), \G\ = 2. So instead of describing the n- 
GHZ state using n generators, we use only one, irrespectively 
of the number of qubits. Thus the rationale behind this ap- 
proach is to find a more convenient description for a certain 
class of stabilizer/graph states. 

In the last section we proved that the generators matrix as- 
sociated to the G-homogeneous state is: 



A k ,n 
On — k.n 



A 



n—k,n 



(45) 



where mn is the mxn zero matrix and A = (ai, .., aj,) T . In- 
tuitively, if we can separate the X- and Z-type generators, we 
can "throw away" the Z part and work only with the spin flips 
(acting on vacuum), since the Z's leave invariant the vacuum 
|0>. 

Consider now some transformations. A Hadamard on the 
i-qubit Hi interchanges Zi <-» Xi, and it is equivalent to ex- 
changing the i-column in the X part with the i column in 
the Z part of the generator matrix j45l >. Therefore applying 
H® n interchanges the X and Z blocks and this is conve- 
nient if n — k < k, since we have a local equivalent state 
described by a smaller group of spin-flips. The next proposi- 
tion establishes the relationship between graph states and G- 
homogeneous states. 

Proposition 4. Let Q be a bipartite graph (i.e., 2-colorable) and 
let V\ , V2 , with I Vl I < I V2 1 , be the sets of vertices forming the 
bipartition (thus all the vertices in Vj have the same colour). 
Then the associated graph state \Q) is locally equivalent to a 
G-homogeneous state, |G) = U\Q), with U = Yiiev 2 ^i- 
Moreover, the group G of spin-flips satisfies |G| = 2^ Vl \. 
Proof. Without loss of generality, we can label the vertices 
such that the first n\ :— |Vi| belong to the V\ partition; define 
also n% := | V2 1 = n — m. Then the generator matrix of the 
stabilizer is 



Or, 



n 



A,. 



0n 2 



(46) 



where the subscripts denote the size of the matrices (for square 
matrices only one index is used); A is the non-zero part of the 
adjacency matrix. Performing a Hadamard on all the qubits 
belonging to V% interchanges the X and Z columns: 



In 



0n 2 



Oni ni ^ n2 
A T 11 



which is of the form (I45> . so we can write 

\g) = n 



(47) 



(48) 
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We can also check that all generators commute, hence 
[X(ai),Z(bj)] = 0, Vi = and j = l,...,n 2 . We 



have (a, , bj 



A,; 



mod 2, as expected (this 



should have been obvious, since applying local Hadamards to 
a generators matrix gives another generators matrix). Since 
the number of X generators is m, we have \G\ = 2l yi L □ 
This result has been independently proved by Chen and Lo 
ll29ll . Moreover, they also proved the reciprocal: any CSS state 
is locally equivalent to a 2-colorable graph state. Therefore we 
conclude that if G is a group of spin flips acting on qubits, the 
corresponding G-homogeneous states are locally equivalent 
to 2-colorable graph states. 

Applying Ea.(ll7>. it is now easy to derive a bound on the 
entropy for bipartite graph states: 



?7 1 

s<io ga iGi = ivii< y 



(49) 



recovering a result obtained in 1 23 ] . 

Examples of 2-colorable graphs include cluster states, trees 
(n-star, n-linear) and 2n-ring graphs. Since they are equiva- 
lent to G-homogeneous states, they have a simpler descrip- 
tion in terms of only Jf-type generators. How widespread 
are 2-colorable graph states? In Ref. 12311 the authors found 
45 (connected) graph states with up to seven vertices which 
are not equivalent (under local unitaries and graph isomor- 
phisms); out of these, 32 are 2-colorable and the other 13 are 
3-colorable. 



C. Qudits 

The formalism of G-homogeneous states can also encom- 
pass multiqudit states. A qudit is a quantum system having a 
rf-dimensional Hilbert space Hd = span{|0), ... ,\d— 1)} = 
C . The generalized Pauli operators for qudits are defined as 



X\k) = \k@l) 
Z\k) = uj k \k) 



(50) 
(51) 



where to — e 2m / d and © is the sum modulo d. They are no 
longer idempotent, since X d = Z = 1; moreover ZX = 
to XZ. It is straightforward to see that for the case of qubits 
(d = 2) we recover the usual definitions and commutations 
relations for X and Z. We can generate any basis vector by 
applying spin flips on |0), \k) — X k \0), hence an arbitrary 

state in Hd has the form = X)fc=o a kX k \0). 

As an example, consider the Hilbert space of n qudits H = 
(C d )®" and let G = (X® n ) = {1, X®"\ (X®™)^ 1 }, 
with X® n — n"=i obviously |G| = d. The correspond- 
ing n-qudit G-homogeneous state is 



|G) =d- 1 / 2 ^(X®") I |0) 



(52) 



i=0 



It is immediate to see that | G) is the maximally entangled state 
of n qudits (it generalizes the n-GHZ state), since d^ = c?b = 
1 for any bipartition (A, B), so the entropy is S — log 2 d as 
expected. 



V. MULTIPARTICLE ENTANGLEMENT 

In this section we show that the compact form of writing 
a G-homogeneous state <[3j is a useful device in calculating 
multiparticle entanglement. The von Neumann entropy char- 
acterizes well the bipartite entanglement of a system in a pure 
state. However, S fails to address the problem of multipar- 
tite entanglement. This becomes apparent even for the simple 
case of three qubits; there are two classes of states with gen- 
uine tripartite entanglement, the W- and GH Z-type states, 
which cannot be distinguished by calculating only the bipar- 
tite entanglement 13011 . An entanglement measure which does 
dis ting uish between the two families of states is the 3-tangle 
t 3 |3l|], since t 3 (GHZ) = 1 and t 3 (W) = 0. More generally, 
Tzi^GHz) > 0, and r 3 (ip w ) = for all the states tp G HZ, ipw 
belonging to the GHZ and the W family, respectively. 

A generalization of the 3-tangle to an even number n of 
qubits is the n-tangle introduced in Ref. 1 32] . For a state 
the n-tangle is defined as: 



\{^\Y^ n \^*)Y 



(53) 



where * means complex conjugation. It has been shown that 
r n is an entanglement monotone and is invariant under local 
unitaries. Note that r n is not defined for n odd. 
Proposition 5. Let G = (X(ai), .., X(afc)) be a group of 
spin flips. Denote by pi := (aj, aj) the parity of the generator 
gi = X(&i) and let p(h) = (a, a) be the parity of an arbitrary 
element h = X(a) 6 G. Then we have the following 

k 

\g\- 1 Yi (- i ) p(,i) = ri( i - = : 1 - p^ g ) < 54 > 

h£G i=l 

where the last equation defines the parity p(G) of the group G. 
The proof is trivial. If all the generators are even, then all the 
group elements have even parity, hence the above sum is 1 and 
p(G) — 0. If only one of the generators is odd (say gi), then 
exactly half of the group elements are odd (those containing 
gi in their expansion) and the other half are even; therefore 
the sum is and p(G) — 1. If more than one generator is odd, 
we can always choose an equivalent set of generators such 
that only one is odd, say gi, e.g., by multiplying all the odd 
generators by g\ (apart from g\ itself); hence this case reduces 
to the previous one. 

We now calculate r„ for G-homogeneous states (and, 
consequently, for 2-colorable graph states, as they are lo- 
cally equivalent). Let = \G\-^ 2 J2 geG g\0) = 

|G|- 1/2 E ae A^( a )|0) = Since F®«X(b) = 

(_l)(b,b)X(b)Y® n = i n (-l)( b - b ) X{b)X® n Z® n and 
Z® n \Q) = |0), we obtain 



(G\Y® n \G) = |G|" 



a,be_4 



X(a)Y® n X(b)\0) 



nor 1 £ 

a.,beA 

(_l)(b,b) ( 



E 

b,b>eA 



= 2 n |G|~ 1 
= i "(l-p(G)) X G(^ n ) 



\X(a)X(b)X® n \0) 
)|X(b')X®"|0) 

(55) 
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where the characteristic function xg{X® u ) = 1 if X® n S G 
and otherwise. Therefore we have the following: 
Corollary. For a G-homogeneous state the n-tangle is 

T n = {l-p{G)) XG {X® n ) (56) 

Hence the n-tangle is 1 iff all the generators of G are even and 

X® n <E G. 



VI. SUMMARY 

Entanglement has been used recently to probe many body 
systems and gain insights of their properties. Quantum phase 
transitions are a notable example in this sense. Several authors 
have shown that a critical behavior of entanglement can signal 
quantum phase transitions ||5[ . 

However, calculating the entanglement entropy for an arbi- 
trary system is often a computationally intractable problem. 
One way to circumvent this problem is to focus on states with 
extra built-in symmetries. In the present article we formal- 
ized this intuition using tools from group theory. We have in- 
vestigated the entanglement entropy for a class of states con- 
structed by acting with the group algebra of a possibly non- 
Abelian group G on a separable reference state |0). The group 
is required to have a bilocal action with respect to a given par- 
tition A, B of the full Hilbert space H = H A ®H B - We 
started first with the G-states, which are constructed as an ar- 
bitrary superposition of group elements acting on |0). These 
states are extremely general, as any state of a Hilbert space can 



be regarded as a G-state with an suitable group G. Moreover, 
we have shown that the ground states of generic Hamiltonians 
are G-states. We have derived an upper bound for the entropy 
provided a separability condition holds for the coefficients. 

A particular class are the G-homogeneous states and in 
this case we generalized our previous results lll9l Eoll . If 
G is a group of spin flips, we show that the associated 
G-homogeneous states are locally equivalent to 2-colorable 
graph states and CSS states. Examples include GHZ states 
and the ground state of the Kitaev model. 

We have shown that we can regard the G-states as the phys- 
ical states in a quantum gauge theory. In this framework all 
the physical states are obtained acting on the G-homogeneous 
state with the commutant algebra of the group algebra of G. 
With some extra assumptions, we can compute the bipartite 
entanglement for all the physical states. We have shown how 
to relate this quantity to the geometric entropy introduced in 
ll33ll . namely the von Neumann entropy relative to a biparti- 
tion obtained by considering a closed surface £ of area a and 
taking as subsystem A all the particles (or degrees of freedom) 
within S. Moreover, the entanglement of the physical states 
obeys the area law, i.e., S{\4> p hy S )) = /(c)- 

Finally, we have shown that this construction can be ex- 
tended in two directions: computing the entanglement entropy 
for qudits in G-homogeneous states and the n-tangle for 2- 
colorable graph states. 

A future challenge will be to find other physical systems 
which can benefit from the group theoretical framework de- 
scribed here. 
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